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ABSTRACT 

The  results  of  coherence  measurements  of  three  free  precession 
proton  magnetometers  are  discussed,,  The  measuring  and  data  proces= 
sing  instrumentation  is  described*  The  noise  associated  "with  the 
instrumentation  has  been  measured  and  presented  in  statistical  form. 
The  results  of  measurements  with  elementary  local  perturbations  are 
included  to  illustrate  the  target  detection  problem*  The  limitations 
of  the  instrument  for  target  detection  are  discussed  and  ir.iporvements 
suggested. 
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1,  Introduction 

The  use  of  magnetometers  for  target  detection  requires  identi- 
fication of  signals  in  the  presence  of  the  nois©  associated  with  the 
earth* s  magnetic  fieldo  The  signal,  in  many  cases,  will  be  smaller 
in  magnitude  than  the  normal  diurnal  variation  of  the  field,  local 
fixed  perturbations  and  geomagnetic  pulsations <>  The  system  noise 
is  of  major  concern  in  a  detection  scheme  using  magnetometers •  This 
paper  presents  the  preliminary  results  of  a  study  of  three  proton 
free  precession  magnetometers  operated  under  non  laboratory  condi- 
tions at  Naval  Postgraduate  School  La  Mesa  Village  sitee 

The  free  precession  magnetometer  consists  of  a  coil  of  wire 
in  which  a  source  of  protons  is  placed „  A  strong  field  is  applied 
to  orientate  the  proton  spins  in  a  direction  approximately  normal 
to  the  earth's  field.  The  polarizing  field  is  then  removed  in  a  man- 
ner such  that  the  spins  cannot  follow  the  fields  The  earth's  field 
then  causes  the  protons  to  precess  about  the  lines  of  magnetic  force 
at  the  Lamour  frequency  given  bys 

uj   =  XPH 

where       ^p    =  gyro  magnetic  ratio 

s  magnetic  moment/angular  momentum        (1) 
and 

W   s  magnitude  of  external  field 

(earth's  field) 

The  net  precessing  magnetic  moment  causes  a  voltage  to  be  induced  in 

a  coil  surrounding  the  sample •  Since  the  processional  frequency  is 

dependent  only  on  the  external  magnetic  field,  the  measurement  of  H 

consists  of  measuring  the  frequency  of  the  induced  emf •  The  inherent 


accuracy  of  the  free  precession  magnetometer  is  determined  by  the 
spectral  line  width  (l)p  which  is  related  to  the  relaxation  time 
T2*  and  is  given  bys 

.  H      _j 

where         --    3.6*  (2) 

H  a  Half  intensity  width 
T«  a  1  sec  (meas) 
s  10°^  GAUSS 
The  gyromagnetic  ratio  has  been  measured  by  the  NBS  to  better  than 
1  part  in  KP„  Full  advantage  of  this  inherent  accuracy  has  not  been 
utilized  because  of  the  problems  associated  with  the  frequency  measure* 
ment.  The  frequency  measurement  is  complicated  by  the  fact  that  only 
a  limited  time  is  available  for  measurement*  The  signal  received 
from  the  magnetometer  is  in  the  few  microvolt  ranges  hence  locally 
generated  electronic  noise  and  transients  from  the  polarizing  opera  <= 
tion  cause  the  signal  to  noise  ratio  to  be  small,,  Standard  counting 
techniques  are  used  for  measuremento  To  get  the  required  accuracy p 
the  frequency  must  be  obtained  from  period  or  multiple  period  average 
measurements o  The  process  is  digital  in  nature  and  lends  itself 
very  nicely  to  computer  processing  of  the  data. 

The  main  disadvantage  is  that  any  reasonable  attempt  to  obtain 
an  analog  signal  proportional  to  the  measurement  is  an  analog  of 
period  rather  than  frequency!  therefore^  the  analog  is  the  inverse 
of  the  quantity  of  interest,.  This  does  not  present  serious  problems 
for  target  detection  since  the  absolute  magnitude  of  field  is  not 
of  interest « 


2.  Instrumentation 
2.1  Measurement 

The  engineering  aspect  of  the  measurement  problem  can  be  stated 
simply  as  follows? 

Given;  Three  transient  type  signals  of  frequency  approximately  two 
kc/s  occuring  at  the  same  time,  with  usable  duration  of  about  3/U 
sec.  The  signals  are  in  the  microvolt  range.  The  recurrance  rate 
is  to  be  2  see*  The  measurements  must  be  conducted  in  the  presense 
of  high  electronic  noise  such  as  relay  transients ,  contact  noise 
and  transients  from  the  collapsing  polarizing  field© 
Finds  The  frequency  of  the  oscillations  accurate  to  1  part  in  10^ 
or  better,  and  suitable  means  of  data  aquisition  and  processing. 

The  funds  available  excluded  the  possibility  of  using  commer- 
cially constructed  instruments  or  componentsc  Also,  modifications 
and  interfacing  would  have  to  be  done  if  components  were  usedo  Hence, 
the  basic  measuring  instrumentation  was  constructed  at  NPGSo  No 
detailed  circuit  design  is  presented^  however  the  circuit  diagrams 
are  included  in  Appendix  IV » 

Assume  the  precession  signal  has  been  amplified  in  a  suitable, 
low  noise,  narrow  band  preamplifier,  and  converted  to  a  signal  hav= 
ing  fast  waveform.  Referring  to  the  simplified  block  diagram  of 
Fig.  1,  the  basic  counter  operation  is  as  follows g  Let  the  slow 
and  fast  counters  consist  of  flip-flop  counting  units9  connected  to 
form  a  binary  counter.  The  gate  control  flip-flop  is  capable  of 
being  triggered  into  either  state.  The  multiplier  converts  the  clock 
frequency  to  a  harmonic  of  the  precision  100  kc/s  oscillator. 


BASIC  COUNTER  UNIT 
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Figure  1 


Let  all  counters  be  reset  and  ready  for  operation,.  The  unknown 
input  signal  is  applied  to  the  gate.,  After  a  delay  of  approximately 
13>0  ms5  an  external  signal  triggers  the  gate  control  to  open  the  sig=» 
nal  gate.  The  long  delay  is  necessary  to  allow  for  damping  the  polari« 
zing  field  and  relay  transients*  The  first  pulse  through  the  signal 
gate  opens  the  clock  gate,,  The  fast  counter  then  counts  the  number 
of  clock  pulses  occuring  during  the  time  interval  determined  by  the 
slow  counter,.  The  slow  counter  is  increased  one  for  each  cycle  of 
unknown  input  signal.  After  the  slow  counter  has  reached  a  predeter~ 
mined  number 9   both  signal  and  clock  gates  are  closed s   and  the  number 
in  the  fast  counter  is  a  value  of  the  tine  interval  of  2n  cycles  of 
unknown  input  signal s   or  can  be  reduced  to  the  average  period  of  the 
signal. 


Let  us  define 

N  -  number  in  fast  counter 

2n  s  number  of  cycles  of  unknown  signal 

tp  s  period  of  clock  ( y^sec.) 

tn  s  period  of  unknown  signal 

Then  we  have 


te/2' 


Substituting  in  (1)  gives 

H   -       *LlLl£  (3) 

and  with  n§10s 

U-     2343400  x/0£4 

M~  /V  (l4) 

The  accuracy  of  such  a  scheme  depends  ong 

a«  The  accuracy  in  determining  the  zero  crossing  of  the  unknown 
signal • 

b»  Propagation  delay  in  the  slow  counter • 

c.  Gating  speed. 

d.  Accuracy  of  the  clock a 

e„     Plus  or  minus  the  las  bit<>  due  to  noncoherent  gating  as 
illustrated  in  Fig.   2. 


EFFECT  OF  NONCOHERENT  GATING 


Sample  Period- 


Figure 2 

The  errors  due  to  determining  the  zero  crossing  are  usually  the  domi- 
nant error  in  using  this  technique!  however  when  averaged  over  102ij. 
cycles  they  become  small  and  the  largest  error  is  the  plus  or  minus 
last  bit,  as  shown.  The  accuracy  in  determining  the  zero  crossing 
in  the  presence  of  noise  (2)  isg 


Ti 


n 
W 


where g 


t  -  maximum  error  in  trigger  time  in  sec, 


;  peak  to  peak  noise  in  volts. 


S  §  slope  of  ideal  signal  at  point  of  trigger  level e 
For  a  sinusoidal  input g 
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For  102h  periods g 

%n-    Q.  3  K  td^seo 

The  plus  or  minus  last  bit  error  if  the  clock  is  100  ke5  is  &-1Q 
sec,  and  the  other  errors  mentioned  are  negligible  compared  to  the 
lOy/f  sec  •  error. 

Consider  the  ^1  bit  limitation  in  detail  for  the  earth's  field 
£0?£0   q    ,  which  is  the  approximate  value  at  Kontereyj,  California* 
lets, 

N  s  OOI3U05I3  <k  1  bit 

1  bit  -^  10^^sec8 
theng 

This  becomes  the  sensitivity  of  the  instrument©  If  the  clock  is 
increased  to  1  mc/sec^  1  bit  s  ^rfc/sec^  andg 

^  H  ■  £  0,11  V'  1  however  other  noise  becomes  appreciable  and 
cannot  be  neglected*  The  circuits  were  designed  for  1  me/sec  opera< 
tion,  but  the  measurements  to  date  were  conducted  at  100  ke/s  and 
300  kc/s  clock  rates.  The  added  complication  is  to  construct  5 
meters  that  will  read  identical^  the  input  circuitry  and  slow  count- 
ers must  have  identical  characteristic s0  This  is  impossible  from 
a  practical  point  of  view,,  The  result  is  that  the  three  counters 
with  an  ideal  signal  input  can  vary  -^  1  bit  in  comparing  their  read' 
ings. 

A   simplified  block  diagram  and  photo  of  the  magnetometer  is 
shown  in  Figs0  3»  h  and  $,     It  must  be  remembered  that  operation 


is  to  be  in  remote  locations  where  extremes  of  temperature,  line 
voltage  fluctuations,  poor  grounds,  etc.  will  be  founds  hence  cir= 
cuits  that  are  sensitive  to  these  variations  must  be  carefully  cora~ 
pensated. 

For  the  experiments  under  consideration,  locating  the  sensing 
coils  approximately  !?00»  apart  would  be  sufficient^  however  for  fur» 
ther  experiments  a  distance  of  greater  than  500 •  might  be  necessary • 
Therefore,  a  provision  to  place  the  preamps  in  the  line  nearer  the 
sensing  coils  for  improved  signal  to  noise  ratio  was  included,,  Bias 
for  the  preamps  is  applied  to  the  signal  line,  from  the  console  to 
eliminate  the  need  for  battery  replacement.  The  relays  that  provide 
the  two  to  three  amps  of  polarizing  current,  and  provide  for  dissapa- 
tion  of  the  stored  energy  when  polarizing  current  is  turned  off,  must 
precede  the  preamp.  Since  this  unit  is  placed  near  the  sense  coil, 
it  must  be  magnetically  clean. 

Referring  to  the  block  diagram  the  components  that  are  common 
to  the  three  units  are  the  clock,  timer,  recorders  and  the  power 
supplies. 
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Figure  5 
Three  proton  magnetometers 

fi7 

and  the  Rb   vapor  magneto- 
meter at  La  Mesa  Village 
Site 
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2*2  Data  Aquisition  and  Processing 

Figure  6  shows  a  block  diagram  of  the  data  collection  and  repro- 
duction system. 

The  scheme  used  for  data  handling  is  certainly  not  optimum  in 
the  sense  of  tape  utilization  and  eat;e  of  operation,  but  an  Ampex 
CF100  was  available  and  required  no  additional  expenditure  for  tape 
recorder.  Also  with  minor  modifications  a  data  recorder  can  be  uti- 
lized at  a  later  date  if  funds  becor.s  available. 

Consider  the  form  of  the  data.  The  sensing  coils  are  polarized 
for  one  second,  at  the  end  of  this  time  the  polarising  current  is 
turned  off,  the  transients  are  dampedj  and  the  signal  is  fed  to  the 
counting  units.  After  102U  cycles  of  the  input  signal  (about  1/2 
sec.)  the  measurement  is  complete.  In  the  fast  counters,  there  are 
three  binary  numbers  to  be  processed.  The  magnitude  of  the  binary 
numbers,  for  H  s  £0950  £  v±%h   clock  frequency  s   100  kc/s,  is  about 
OOI3I4IU60  or  16  bits.  With  the  clock  increased  to  300  kc/s,  the 
number  is  001*25 260g  or  18  bits.  Hence,  at  the  end  of  the  counting 
interval,  there  are  US  bits  to  be  recorded  in  a  format  for  computer- 
input.  The  Ampex  CP100  is  a  "portable"  instrumentation  recorder 
with  lU  parallel  inputs  and  capable  of  speeds  from  1  7/8  ips  to  60 
ips.  All  of  the  record  and  reproduce  amplifiers  that  were  on  hand 
were  not  capable  of  recording  D.C.  level sj  therefore,  the  information 
was  put  on  tape  in  the  form  of  pulses.  Also,  for  reasonable  speed 
in  reproduction,  the  data  was  recorded  at  1  7/6  ips  and  reproduced 
at  60  ips,  A  record  of  approximately  6*$   hours  can  be  reproduced  in 
10  minutes.  The  machine  used  to  read  the  raw  data  is  a  CDC160  com- 

12 


DATA  AQUISITION  AND  PROCESSING 
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I 

Graph  Plotter 

Print 

puter  capable  of  reading  words  of  up  to  12  parallel  bits.  The  mag- 
netometer words  are  broken  into  two,  12  bit  words  and  time  multiplexed 
on  tape  to  satisfy  the  computer  input  conditions6  Hence ,  for  each 
sample,  there  are  six  12  bit  words  to  be  processed.  Even  with  the 
time  speed  up  in  reproduction  the  CDC  160  is  operating  below  its 
capabilities,,  In  an  effort  to  utilize  the  160  time,  the  data  was 
put  on  tape,  so  that  when  it  occured,  the  160  would  read  at  its  maxi- 
mum  rate  and  use  the  Mno  data  time"  to  continue  with  preliminary 
processing*  Channels  13  and  Ik  on  the  tape  are  input  ready  and  dis~ 
connect  signals  used  to  control  the  l60„  And  of  course,  for  reason- 
able machine  time  efficiency,  there  is  on  hand  enough  tape  to  allow 
continuous  recording  for  one  week. 

The  Ampex  CP1C0  was  not  designed  specifically  for  this  type  of 
operationj  hence  the  need  for  an  interface  to  mate  it  with  the  l60o 
The  interface  must  provide  the  following  functions? 

a.  Convert  the  pulse  type  data  to  D.C,  level  signals. 

b.  Put  the  data  on  line  when  a  request  is  supplied  from  the 
computer, 

c.  Effectively  disconnect  itself  from  the  lines  when  other 
equipment  is  called© 

d.  Provide  control  information,  i„e.  input  ready  and  discon<= 
nect  for  the  160„ 

e.  Effectively  isolate  the  computer  from  spurious  transients 
from  the  tape  machine „  These  normally  occur  during  start  «  stop 
operations,  but  the  recorder  is  somewhat  susceptible  to  power  line 
transients  and  they  sometimes  appear  in  its  outputo 


1U 


Figures  8,  9  and  10  shows  the  block  diagram  of  the  interface 
provided,  Figo  11  is  a  photo  of  the  preliminary  processing  setup. 
The  cards  were  assembled  and  mounted  as  a  component  on  a  laboratory 
cart,  the  tape  unit  is  placed  on  the  top  of  the  cart  and  computer 
cable  connectors  provided.  This  eliminates  the  need  for  a  permanent 
rack  in  the  computer  laboratory • 
The  preliminary  processing  consists  of  the  following g 

a,  Reading  the  raw  data, 

b,  Checking  the  data  and  keeping  a  record  of  the  errors, 

c,  Writing  the  storage  tape. 

d,  Froviding  identifying  blocks  on  the  storage  tape. 

e,  Searching  the  storage  tape  for  any  particular  record  or 
the  end  of  the  data  so  that  one  storage  tape  may  be  used  for  many 
week's  records, 

f ,  Direct  display  of  the  raw  data. 

The  data  storage  format  used  is  illustrated  in  Fig.  7* 

STORAGE  TAPE  FORMAT 
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22  BLOCKS 
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Figure  7 


With  the  time  base  ussd  there  are  nearly  380l|G  twelve  bit  words 
per  6.5  hours  record.  The  6.5  hours  record  is  one  3600  roll  of  in- 
strumentation tape  at  1  7/8  ips.  Since  the  operation  is  stopped  to 
change  tape,  it  is  convenient  to  put  the  identifying  block  at  the 
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end  of  this  record.  Also,  this  number  of  data  words  can  be  proces- 
sed in  the  order  of  minutes  by  the  CDC160U  so  that  if  the  computer 
is  needed  for  other  jobs|  it  can  be  stopped  at  the  identifying  blocks 
and  return  later  to  the  same  point  with  the  proper  index  setting. 
The  block  length  was  set  somewhat  arbitrarily,,  The  memory  storage 
limit  of  the  160,  write  time  of  the  163  tape  unit  and  suitable  re= 
cords  for  handling  in  the  160U,  are  the  factors  to  be  considered  when 
chosing  the  block  length.  The  instrumentation  recorder  cannot  con- 
veniently be  stopped  when  writing  the  storage  tape 5  hence  one  data 
point  is  lost  when  transferring  data.  This  represents  22  points 
per  6.5  hours  record  or  less  than  a  minute  in  real  time.  The  block 
length  was  chosen  to  be  3li56io  twelve  bit  words  per  block 9  the  fig- 
ure being  compatable  with  the  above  conditions,, 

The  identifying  record  consists  of  writing  the  months  day  and 
year  in  flex  code.  The  remaining  cells  of  the  block  are  loaded  with 
zeros.  It  is  then  easy  to  separate  the  data  from  the  identifying 
blocks  since  none  of  the  data  v<ords  will  be  all  zeros* 

The  fact  that  the  data  had  to  be  reconstructed  in  the  loOU  was 
also  considered  in  chosing  the  block  length.  The  information  is 
buffered  into  the  l60l|  in  U8  bib  words.  The  blocks  consist  of  an 
even  number  of  data  points,  i.e.  six  12  bit  or  two  I4.8  bit  words t 
After  a  block  is  read  into  the  machine  and  loaded  into  the  A  register, 
the  magnetometer  readings  can  be  stored  sequentially  in  memory  with 
a  loop  of  a  few  machine  language  instructions.  The  data  analysis 
is  then  done  using  Fortran,  allowing  easy  and  versatile  manipulation 
of  the  data.  The  160  and  160U  input  programs  are  included  in  Appen= 
dices  II  and  III. 
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Reassembly  is  illustrated  below. 


160U  WORDS 


Mag  1  (2U  bits) 

Mag  2  (2U  bits)   | 

Mag  3  (2li  bits) 

Mag  1  (2U  bits) 

Mag  2  (2U  bits) 

Mag  3  (21*  bits) 

Briefly,  the  operation  is  as  follows  %     Magnetometer  readings 
are  recorded  for  long  periods  of  time  at  a  remote  site  on  a  pro table 
magnetic  tape  unito  The  tape  unit  is  then  returned  to  the  process- 
ing site,  and  using  a  time  speed  up  of  X  32,  the  intermediate  pro= 
cessing  is  performed  and  the  data  stored. 

The  storage  tape  along  with  suitable  programs  is  then  placed 
on  file  for  processing  by  the  computer  facility  at  their  convenience. 
The  tape  unit  is  then  returned  to  magnetometer  site  for  more  records. 
The  new  data  obtained  is  added  to  the  storage  tape  so  as  not  to  de» 
stroy  the  existing  data. 
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3»  Data  Analysis 

The  experimental  measurement  program  was  set  up  to  study  the 
coherence  between  the  magnetometers^,  i,e<>  the  comparison  of  simulta- 
neous samples  of  the  field,  with  the  sensors  located  at  different 
spatial  locations.  The  coherence  was  to  be  examined  with  various 
local  perturbations  in  the  field.  With  the  target  detection  objective 
in  mind,  the  least  complicated  method  of  obtaining  the  target  signal 
will  result  in  the  best  instrumentation  for  field  use.  The  problem 
is  somewhat  analogous  to  the  radar  problem  where  decisions  must  be 
made  in  the  presence  of  noise,  with  associated  probabilities  of  error 
in  detection. 

There  are  several  types  of  geomagnetic  noise  that  occur*  The 
most  dominant  is  the  diurnal  variations.  Superimposed  on  this  are 
smaller  variations  known  as  raicropulsations  and  noise  caused  by  lo- 
cal disturbances.  Sudden  eruptions  on  the  sun  sometimes  cause  large 
fluctuations  in  the  field.  Several  studies  of  the  various  noises 
have  been  made  (3),  (U)  an^  ($)•     For  the  proton  precession  magneto- 
meter, the  instrument  produces  noise  of  major  contribution  to  the 
spectra  of  small  geomagnetic  pulsations.  For  target  detection^  the 
signal  of  interest  will  most  likely  be  of  the  same  order  of  magni- 
tude as  the  small  micropulsations  and  local  distrubanees.  For  geo=> 
physical  prospecting  and  mapping,  the  signals  will  be  much  larger,. 

As  with  any  raw  data,  there  will  be  some  points  that  are  obvi- 
ously in  error.  The  cause  of  errors  is  usually  uncertain^  but  some 
provision  must  be  made  in  the  processing  to  eliminate  these  errors 
so  as  not  to  distort  the  data.  Care  must  be  taken  when  performing 
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these  operations  so  the  target  signals  are  not  masked.  Records  of 
the  errors  should  be  kept  in  order  that  decisions  about  the  useful- 
ness of  the  set  of  data  can  be  made. 

The  problem  of  coherence  between  the  spatial  samples  requires 
that  the  instrument  noise  be  determined^  with  an  ideal  signal^  then 
measurements  be  conducted  with  no  local  perturbations  in  the  fields 
and  the  samples  compared,,  The  simplest  comparison  is  to  plot  the 
difference  signal  vs.  time.  With  no  local  perturbation^  this  means 
essentially  a  plot  of  the  system  noise  vs0  timej  hence  the  need  for 
filtering  is  obvious.  To  use  this  method  for  target  detection  will 
require  serious  consideration  of  the  filter  time  constants* 
The  types  of  filtering 

1,  A  clipping  filter  (on  the  raw  data)  that  clipped  at  4"  200  yc 
The  number  of  points  outside  of  this  range  are  recorded *  This  range 
was  about  twice  the  expected  diurnal  variation*  The  filter  removes 
the  errors  from  lost  high  order  bits  and  other  instrument  sources. 
The  typical  error  rate  was  about  1/50000. 

2«  A  simple  three  point  average  for  display  of  the  raw  data  on 
the  DD  65  display  unit  was  used.  This  gives  an  operator  a  fast  vis= 
ual  display  of  the  raw  data.  Decisions  can  then  be  made  about  the 
quality  of  the  data. 

3.  Eighteen  point  averaging,  i.e0  36  seconds  real  time5  gave 
partially  smoothed  6.5  hours  graph  records, 

U.  A  simple  single  section  low  pass  filter  was  used  of  the 
form  l/(s  °  a).  The  time  constant  was  set  at  100  sec  and  20  min, 
in  real  time.  The  implementation  was  done  using  Z  transform  (6). 
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The  filter  design  was  done  by  Professor  Harold  Ao  Titus. 


We  have 


*in(S) 


-  a 


which  has  a  Z  transform 

X„C2)     ^jQ 

In  block  diagram  form  this  is  seen  in  £ig<,    21 


Mli)J 


XgU 


i 


£ 


■*- 


Ko(k) 


Figure  29 
In  terms  of  the  discrete  time  samples  we  have 
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a.  /vt^p 
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T  ~  sample  period 

a  s  time  constant 
The  function  was  trunkated  at  n  §  20  for  values  of  k  ^  20„ 
The  contribution  to  the  filtered  value  from  the  terms  multiplied 
by  e"anT  for  n  >  20  is  negligible. 
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The  coherency  is  a  statistical  parameter s  and  must  be  repre~ 
sented  as  such. 
Lets 

/\  1  =  Hl  -  H2 
A  2  =  %  -  H3 

A  3  ■  H2  -  »3 

The  differences  of  the  unfiltered  values  were  grouped  in  a  frequency 
distribution  which  approximates  the  continuous  density  function  for 
large  sample  size*  Since  the  process  is  digitals  the  allowable  dif- 
ferences, in  bits  ares  0,  ^  1,  |  2S   £  °o°.  Therefores  the  frequency 
distribution  is  not  a  smooth  curve.  Since  the  instrument  noise  is 
randomly  distributed ^  the  envelope  of  the  distribution  is  a  normal 
density  function,,  The  mean  was  removed  before  plotting  9   and  the 
statistic  representing  the  standard  deviation  was  computed  according 
tos 


iJc^n 


N 

1" 

where  %  N  a  the  number  of  samples  in  a  6«5  hour  record, 


Comparison  of  the  density  function  curve  with  an  ideal  signal 
input  to  that  of  magnetometer  input  gives  a  measure  of  the  noise 
produced  by  the  device. 

To  further  define  the  coherence  s   the  auto  and  cross  correla=> 
tion  functions  were  computed  and  the  coherency  computed  as  defined: 
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■where: 
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This  discrete  computational  form  is 
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where i 


N  g  number  of  data  values 
n  -  number  of  lags 

Since  the  functions  are  from  finite  time  series  nrnAX  can  only  be  made 
about  Ool  N»     Tha  6e5>  hour  records  vrill  vary  in  form  depending  upon 
the  time  of  day  the  records  were  made.     Therefore 9   the  correlation 


functions  will  not  be  the  same  shape  for  each  record;  however  the 
coherence  function  as  defined  will  be  the  sames  and  for  perfectly 
coherent  records  is  unity „  The  programs  are  included  in  Appendix 
III. 

It  should  be  pointed  out  that  when  perturbations  are  placed  in 
the  field,  the  basic  assumptions  regarding  correlation;  i0e0  a  sta- 
tionary random  process,  are  no  longer  applicable 0 

Experimental  data  was  taken  to  determine  typical  signatures  for 
elementary  perturbations,  Fig»  12  shows  the  geometry  of  the  detector 
centered  coordinate  system,, 
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The  dipole  moment  was  orientated  parallel  to  the  earthvs  field  and 
its  motion  was  in  the  x  direction  such  that  dx/dt  g  constant 0     The 
approximation  r  »  1  is  valid » 
At  the  detectors 


And  j 


H  3    \T  Ho2  °  Hp2 


HpsKM 

13 


where s 

M  s  dipole  moment 

r  s  function  of  time  and  the  geometry 
K  s  constant 

No  attempt  was  made  at  this  time  to  compare  the  measured  values  with 
the  theoretical  values s   because  the  terrain  does  not  lend  itself 
keeping  L$  d,  and  v  constant.  The  results  are  representative  of 
typical  signatures  that  vnsuld  be  found  in  practice e 


28 


U.  Results 

1^.1  Instruments,  ticns 

An  ideal  signal  was  substituted  for  the  procession  signal  and 
set  of  records  were  made,  as  shown  in  Figs,  13  through  2S>„  The 
graph  scales  were  left  unchanged  because  it  was  desired  to  show  the 
instrument  errors  on  the  same  scale  as  the  magnetic  field  scales* 
The  magnetic  field  vs,  time  graphs  were  not  included  because  the  in- 
put signal  was  an  audio  oscillator  and  only  a  small  drift  puts  it 
off  the  graph  scale*  The  clock  was  300  kc/s  which  makes  1  bit  a 

0.36  y  . 

Figure  28  shows  the  precesion  signal  from  the  various  sensors. 
Figs,  26  and  27  are  typical  plots  made  from  the  analog  section  of  the 
instrument.  For  comparison,  the  Kb"'  magnetometer  output  was  adjust^ 
ed  to  the  same  scale  and  plotted  on  the  same  graph.  The  Kb"'  magneto- 
meter output  is  an  analog  of  frequency.,  while  the  proton  magnetometer 
the  analog  is  of  period^  hence  the  reason  for  the  opposite  drift  on 
the  chart.  Figure  27  is  an  analog  record  of  the  differences,,  The 
filter  was  a  single  section  RC  filter  with  time  constant  g  10  sec. 
On  the  multiple  curve  plots  showing  coherence,  the  curves  are  identic 
cal  and  cannot  be  distinguished  from  each  other.  These  curves  can 
be  compared  with  the  curves  obtained  from  the  sensor  with  no  local 
perturbations  in  the  field,  for  the  contribution  of  the  measuring 
instrument  noise  to  the  total  noise. 

Difficulties  with  both  instrumentation  and  computer  tape  have 
been  encountered.  The  use  of  new  computer  tape  for  data  storage  and 
scratch  tapes,  and  frequent  head  cleaning  of  the  instrumentation 
recorder  reduced  the  difficulties  to  a  minimum,, 
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Differences  in  magnetic  intensity  vs.  time 

for  an  ideal  signal  input. 

DELI  a  Hi  -  H2 

DEL2  3  Hi  -  H3 

DEL3  =  H2  -  H3 

Filter  3  18  point  average 
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Figure  13 
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Differences  in  magnetic  intensity  vs. 

time  for  an  ideal  signal  input. 

DEL  1  3  Hi  -  H2 

DEL  2  =  Hi  -  H3 

DEL  3  s  H2  -  H3 

Filter  a  Z  transform  low  pass  with 

time  constant  of  20  minute s* 
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Figure  ll; 


X-SCfiLE  -  L^L+m  UNITS/INCH. 
Y-SCfiLE  -  2.eG£fc£&JJUT3.--INCH. 

PERSON  FILTER 

TLTEREDDIFFERENCE  IN  MflG  FIELD  T  IN  HRS  H  GAMMA 


31 


Normalized  density  function  for  the 
difference  Hi  -  H2  -  mean  (Hi  -  H2)> 
with  an  ideal  signal  input. 
Mean  3  0.061  gamma 
Standard  deviation  =  0,21 
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Figure  15> 
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Normalized  density  function  for  the 
difference  H]_  -  H3  -  mean  (Hi  -  H3), 
with  an  ideal  signal  input. 
Mean  a  0.056  gamma 
Standard  deviation  -  0.72 
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Figure  16 
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Normalized  density  function  for  the 
difference  Ho  -  H3  -  mean  (H2  -  H3), 
with  an  ideal  signal  input. 
Mean  ■  -0.005S  gamma 
Standard  deviation  -  0.72 
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Figure  17 
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Auto  correlation  for  an  ideal 

signal  input. 

Record  length  =  6.5  hours 

Hi  r  Correlation  (Hi  -  mean  Hi) 

H2  =  Correlation  (j^  -  mean  H2) 

H3  s  Correlation  fa  2  -  mean  H3) 

Lag  =  A  t 
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Cross  correlation  for  an  ideal 

signal  input. 

Record  length  =  6.5  hours 
Hi 2  »  Correlation  (Hi  -  mean  Hi  and  H2  -  mean  H2) 
H13  =  Correlation  (Hi  -  mean  Hi  and  H3  -  mean  H3) 
H23  3  Correlation  (Ho  -  mean  H?  and  H3  -  mean  H*0 
Lag  s£  t  *  *' 
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Figure  19 
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Coherence  for  an  ideal  signal  input 
Record  length  :  6,$  hours 
Hi 2  =  Coherence  of  Hi  and  H2 
H13  =  Coherence  of  Hi  and  H3 
H23  s  Coherence  of  H2  and  H3 
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Figure  20 
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Differences  in  magnetic  intensity 
for  an  ideal  signal  input 
Record  length  =  1  hour 
Filter  *   2  point  average 
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Figure  21 
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Difference  in  magnetic  intensity 
for  an  ideal  signal  input 
Record  length  =  1  hour 
Filter  =  Z  transform  low  pass 
with  time  constant  of  100  seconds. 
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Figure  22 
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Auto  correlation  for  an  ideal 

signal  input. 

Record  length  -  1  hour 

Hi  3  Correlation  (Hi  -  mean  Hi) 

H2  s  Correlation  (H2  -  mean  H2) 

H3  3  Correlation  (H2  -  mean  H3) 

Lag  s^t 
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Figure  23 
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Cross  correlation  for  an  ideal  signal  input. 
Record  length  2  1  hour 

Hi 2  *  Correlation  (Hi  -  mean  Hi  and  H2  -  mean  H2) 
Hi 3  «  Correlation  (Hi  -  mean  Hi  and  H3  -  mean  H3) 
H23  =  Correlation  (H2  -  mean  H2  and  H3  -  mean  H3) 
Lag  a  A  t 
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CROSCORREUTTION  FUNCTION  Y  IN  PRODUCTS  X  IN  LAGS 
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Coherence  for  an  ideal  signal  input 

Record  length  =  1  hour 

H12  =  Coherence  of  tt±   and  Hg 

H13  3  coherence  of  H]_  and  H3 

H23  s  Coherence  of  H2  and  H3 
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Figure  25 
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TYPICAL  ANALOG  PLOT  OF 
MAGNETIC  INTENSITY  VS.  TIME 
FILTER  m   RC  WITH  TIME  CONSTANT  OF  6  SECONDS 
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TYPICAL  ANALOG  PLOT  OF 

THE  DIFFERENCES  VS.  TIME 

FILTER  m  RC  WITH  TIME  CONSTANT  OF  6  SECONDS 
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FIGURE  27 
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It* 2  Magnetic  field  measurements 

Only  representative  results  of  the  magnetic  field  measurements 
are  presented^  since  the  measurements  are  not  complete  at  this  time. 
Samples  of  plots  made  with  the  clock  at  100  kc/s  are  included  to 
illustrate  the  effect  of  the  clock  frequency  on  the  results.  For 
the  associated  density  function  of  the  differences $  the  allowable 
values  are  0  *>  1.1,  2  2o2j  *-°  °  <>  in  gamma 9  with  the  clock  at  100 
kc/s.  The  expanded  difference  vs.  time  shows  clearly  the  *  last 
bit  limitation.  Comparing  the  density  functions  with  the  plots  for 
the  clock  at  300  kcg  shows  that  most  of  the  noise  associated  with 
the  measurement  is  due  to  the  last  bit  and  can  be  reduced  by  increas- 
ing the  clock  frequency.  The  perturbation  experiments  were  conducted 
with  the  clock  at  300  ke/s9  The  results  of  the  perturoataon  experi- 
ments show  clearly  the  effect  of  the  digital  filtering 9   and  the  im- 
portance of  tame  constant  consideration  for  target  detection.  The 
oscillator  is  to  be  run  1  mc/s  for  further  experiments*  Also*  many 
of  the  errors  that  appear  in  the  data  were  the  direct  result  of  the 
instrumentation  tape.  New  tape  has  been  purchased  and  should  elim- 
inate most  of  these  errors^,  in  future  records. 

As  pointed  out  previously 9   the  basic  assumptions  regarding  cor= 
relation  are  not  valid^  when  the  noise  is  not  stationary;  however s 
since  correlation  detectors  can  be  instrumented  and  may  prove  useful 
in  the  detection  problem.  Some  of  the  correlation  function  plots  are 
included  to  illustrate  the  effect  of  local  perturbations  on  the  func= 
tions. 
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Jj>«  Conclusions 

For  target  detection^,  where  the  signals  under  consideration,  are 
in  the  \  gamma 5  or  less  range^  it  will  be  difficult  to  detect  bhe 
target  without  resorting  to  sophisticated  signal  processing*  In 
this  range  of  signal 3 s   consideration  should  be  given  to  other  instru- 
ments! such  as  the  Rubidium  optically  pumped  magnetometer 9  whose 
theoretical  sensitivity  is  0o001  gamma 9   especially  where  the  detector 
is  stationary  #  and  orientation  is  not  a  problems  It  can  be  demonstrated 
that  when  using  the  same  measurement  technique  as  for  the  proton 
magnetometer,  the  instrument  and  processing  noise  will  be  an  order  of 
magnitude  less  than  the  theoretical  limits  imposed  by  bhe  sensiti- 
vity  of  the  instrument. 

For  signals  in  the  J  to  if  gamma  range s   detection  could  be  ac- 
complished with  analog  type  correlation  detectors..  These  could  be 
implemented  with  operational  amplifiers. 

Signals  above  this  range  could  be  detected  with  the  simple  dif- 
ferencing technique  $  and  the  associated  R0C0  filters,, 

Increasing  the  clock  frequency  and  averaging  the  period  over 
more  than  102U  cycles  would  reduce  the  measurement  noise  to  a  smaller 
part  of  the  total  noise.,  The  target  signals  could  then  be  reduc  6 
by  a  factor  of  three ^  with  no  change  in  the  instrumentation* 

Full  utilization  of  the  inherently  narrow  line  width  of  proton 
magnetometer  cannot  be  achieved  using  the  standard  measuring  tech- 
nique* Other  techniques  have  been  suggested  by  Professor  Carl  i„ 
Kennekenj  such  as  utilising  a  rapid  coherent  polarization  in  a  phase = 
lock  loop  arrangement*  Further  development  work  needs  to  be  don-3 
to  explore  these  possibilitieso 
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with  20  gamma  step  function  applied  to  H3. 

Record  length  =  6*$   hours 

Clock  =  300  kc/s 

Filter  =  Z  transform  low  pass  with 

time  constant  of  20  minutes 

DEL  1  =  Hi  -  H2 

DEL  2  =  Hi  -  H3 

DEL  3  =  H2  -  H3 


■U3 


&Ji 


£^ 


^ 


J3EL1 


X  SCALE  =  1.00E-+00  UNITS/INCH 
Y  SCALE  -  5.00E400  UNITS/INCH 

FILTERED  DIFFERENCE  IN  MAG  FIELD  T  IN  HOURS  H  IN  GAMMA 
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Earth's  magnetic  field  vs.   time 

Filter  r  13  point  average 

Clock  =  300  kc/s 

Record  length  a  6»$  hours 

Hi  =  H|_  -  mean  Hi 

H2  =  H2  -  mean  H^ 

H3  s  H3  -  mean  Hi 

Perturbations : 

i:  5  gamma  and  +■  10  gamma  applied 

at  constant  velocity  to  H3.   -  2 

gamma,  -  5  gamma,  -  5>  gamma  applied 

at  constant  velocity  to  Ho,  H«  and 

Hi  during  first  hour,     A  5  gamma 

step  applied  to  Ho  at  t  =  1  hour  and 

removed  at  t  -  1  j/U  hours. 


X-SCfLE  -  1.8GE+Q0  UNITS/INCH. 
Y-SWLE  -  L8GE+81  UNITS/INCH. 


ANDERSON  BOX  263 
EARTHS  MAGNETIC  FIELD  US  TIME 


T  IN  HRS  H  GflMMf 
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Earth's  magnetic  field  vs.  time 

Filter  a  Z  transform  low  pass  with 

time  constant  of  20  minutes 

Clock  =  300  kc/s 

Record  length  ;  6,5  hours 

Hi  -  H^  -  mean  & 

H2  ;  H„  -  mean  HJ 

H3  3  Ho  -  mean  H^_ 

Perturbations : 

±  $  gainma  and  +  10  gamma  applied 

at  constant  velocity  to  H3.  -  2 

gamma,  -  5>  gamma,  -  5  gamma  applied 

at  constant  velocity  to  Ho,  Ho  and 

Ht   during  first  hour.     A  5  gamma 

step  applied  to  Ho  at  t  =  1  hour  and 

removed  at  t  :  1  j/h  hours. 


X-3CflLE  -  JL88E+08  UNITVINCH. 
Y -SCALE  «  1.0QE+01  LiNITS/'INCH. 


ANDERSON  FILTER 
EARTHS  MAGNETIC  FIELD  US  TIME 


T  IN  HRS  H  GflMM 
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Earth's  magnetic  field  vs.  time 
Filter  -  2  point  average 
Clock  r  300  kc/s 
Record  length  =  1  hour 

Hl  =  %  "  mean  Hl 

H2  z   H«  -  mean  H]_ 

H3  =  Ho  -  mean  H^ 

Per tur cations: 

♦  5  gamma  and  ±  10  gamma  applied 

at  constant  velocity  to  Ho«  -2 

gamma,  -5>  gamma,  -5  gamma  applied  at 

constant  velocity  to  Ho,  H2  and  H]_ 

during  first  hour.  A  5  gamma  step 

applied  to  Ho  at  t  3  1  hour  and 

removed  at  t  =  1  3/4  hours. 


lrvP 


X  SCALE  =  1.00E+00  UNITS/INCH 
Y  SCALE  -  1.00E401  UNITS/INCH 

EARTHS  MAGNETIC  FIELD  VS  TIME  T  IN  MINUTES  H  IN  GAMMA 
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Earth's  magnetic  field  vs.  time 
Filter  r  2  transform  low  pass  with 
time  constant  of  100  seconds,  ' 
Clock  a  300  kc/s 
Record  length  -  1  hour 
H^  r  H-j_  -  mean  H-j 
H2  r  Hp  -  mean  H-j 
H3  r  H3  -  mean  Hj_ 
Perturbations: 

z   5  gamma  and  =  10  gamma  applied 
at  const-ant  volocity  to  Ho.  -  2 
gamma,  -  5  gamma,  -  $  gamma  applied 
at  constant  velocity  to  H3,  H9  and 
#L  during  first  hour.  A  5  gamma 
step  applied  to  Ho  at  t  -  1  hour 
and  removed  at  t  i   1  3A"*hours. 


X  SCALE  =  1.00E-00  UNITS/INCH 
Y  SCALE  =  1.00E-01  UNITS/INCH 

EARTHS  MAGNETIC  FIELD  VS  TIME  T  IN  MINUTES  H  IN  GAMMA 
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HC 


Auto  correlation  for  earth's  magnetic 

field  vs.  time. 

Record  length  ■  1  hour 

Clock  =  300  kc/s 

H-j_  =  Correlation  (H^  - 


hi 


mean  Hi) 
(Hj  -  mean  H2) 
(H3  -  mean  H3) 


H2  r  Correlation 

Ho  -  Correlation 

Lag  ::  A  t 

Perturbations': 

±  5   gamma  and  i  10  gamma  applied  at 

constant  velocity  to  H3.  -2  gamma, 

-$  gamma,  -5  gamma  applied  at  constant 

velocity  to  H3,  H2  and  ^  during  first 

hour.  A  5  gamma  step  applied  to  K3  at 

t  -  1  hour  and  removed  at  t  =  1  3/u  hours. 


H3 


081 


002 


803 


084 


055 


CC- 


007 


X  "SCALE  -  LOOE+01  UNITS/INCH. 
Y-SC3LE  -  1.0OE+80  UNITS/INCH 


IDERSON  BOX  263 

ITOCORRELflTION  FUNCTION  Y  IN  PRODUCTS  X  IN  I  AGS 


APPEITOIX  I 


H23 

mi 


Cross  correlation  for  earth's  magnetic 

field  vs.   time. 

Record  length  -  1  hour 

Clock  -   300  kc7s 

H-j_p  z   Correlation  (H-,  -mean  K]_  and  HL  -mean  H2) 

H13  r  Correlation  (Hf  -mean  H^  and  Ho  -mean  H3) 

H23  r  Correlation  (H2  -mean  H2  and  H3  -mean  H3) 

Lag  :^t 

Perturbations : 

±  5   gamma  and  ±   10  gamma  applied  at 

constant  velocity  to  Ko«  -2  gamma, 

-5  gamma,  -$   gamma  applied  at  constant 

velocity  to  H3,  H?  and  H]_  during  first 

hour.  A  $   gamma  step  applied  to  H3  at 

t  s  1  hour  and  removed  at  t  =  1  3/u  hours. 


H13 


631 

002 

883 

mi 

835 

GOB 

C37 


X -SCALE  -  L88E+81  UNITS/INCH. 
Y-3CR.E  -  1.88E+88  UNITS/INCH. 


PERSON  BOX  263 

^CORRELATION  FUNCTION  Y  IN  PRODUCTS  X  IN  LfiGS 
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Coherence  for  earth's 

vs •  time . 

Record  length  =  1  hour 

Clock  =  300  kc/s 

Ht  2  r  Coherence  of  H]_  and  IU 

H^3  =  Coherence  of  Hj_  and  H3 

H23  -  Coherence  of  H2  and  H3 

Lag  =  A  t 

Perturbations: 

•4-  5  gamma  and  ±  10  gamma  applied  at 

constant  velocity  to  H3.     -2  gamma, 

-S>  gamma,  -5  gamma  applied  at  constant 

velocity  to  H3,  H0  and  H^  during  first 

•hour.     A  5  gamma  step  applied  to  Ho  at 

t  =  1  hour  and  removed  at  t  =  1  3/u  hours, 


881 


C32 


£33 


69^ 


085 


006 


00" 


X-SCflLE  -  1.8aE+01  IMTS/IMCH. 
Y-SCflLE  -  Q.88E-81  UNITS/IN^ 


ANDERSON  BOX  263 

COHERENCE  FUNCTION  Y  IN  C0H  X  IN  LAGS 
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OPERATOR  INSTRUCTIONS 
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160  PROGRAMS  AND  OPERATOR  INSTRUCTIONS 


Starting,  Addresses 

0000 
0002 
0004 
0006 
0010 
0100 


Function 

Rewind 

Read 

Search 

Write 

Write  identification  Block 

Main  Input  Program 


ERR01 
ERR02 
ERR04 
ERR03 


Cause 

773  of  current  input  in  error 
Memory  storage  not  complete 
End  of  storage  tape 
163  not  ready 


Addresses  of 
Important  Counters 

0057 
0045 
0046 


Information 

Number  of  blocks  written  on 
storage  tape 

Current  storage  location  in 
memory 

Number  of  zero  words  in  mem- 
ory 


Operator  Instructions 

1.  Load  program  and  make  the  necessary  cable  connections, 

2.  Set  identifying  records  in  cells  30q  to  37g° 

3.  Write  identifying  block  (run  from  P  s  010) .   Stop  at 
P  -  100. 

4.  Set  cell  100  -  2200 

57  3  0 

5.  Start  Ampex  C.P.  100  and  run  from  P  -  100.   Stop  when 
tape  is  finished  or  ERR  stop. 


II-l 
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INPUT  PROGRAM 


5f  AST 


Aa  Oy 

/ > 

- 

OUHP 

S8  ft 

ftl 

>s>a 

v^MD 


APp« 

s/ 000 

A 

CtJVRl 

5« 

-  0 

/> 

•w 

s  Mo 

CELL 

=  0 

v 

£AU 

rwpvr 

^ 

i 

••  <■ 

YES 


IWOR 


II-2 


0100 
0101 
0102 

0105 
0104 

0105 
0106 
0107 
0110 
0111 
0112 

0113 
0114 

0115 
0116 
0117 
0120 
0121 
0122 
0123 
0124 

0125 
0126 
0127 
0130 
0131 
0132 
0133 
0134 

0135 
0136 

0137 
0140 
0141 
0142 
0143 
0144 
0145 
0146 
0147 
0150 
0151 
0152 

0153 
0154 

0155 
0156 


APPENDIX  II 

2200 

LDF00 

1000 

1000 

4045 

STD45  (ADDR) 

0400 

LDNOO 

4046 

STD46  (CNTR1) 

0400 

LDNOO 

4150 

ST  1 50 

4151 

STI51 

4152 

ST  1 52 

4153 

STI53 

4154 

ST  1 54 

4155 

ST  1 55 

7500 

SXFOO 

0500 

0500 

7203 

INP03 

0620 

0620 

6102 

NZF02 

0600 

0600 

1600 

LSFOO 

0607 

ADN07 

6115 

NZF15 

2150 

LDI50 

6013 

ZJF13 

2151 

LDI51 

6011 

ZJF11 

2152 

LDI52 

6007 

ZJF07 

2153 

LDI53 

6005 

ZJF05 

2154 

LDI54 

6003 

ZJF03 

2155 

LDI55 

6105 

NZF05 

5446 

A0D46 

0777 

SBN77 

6536 

NZB36 

0001 

ERR01 

2050 

LBD50 

4047 

STD47 

2147 

LDI47 

4145 

ST  1 45 

5445 

A0D45 

5447 

A0D47 

3456 

SBD56 

6505 

NZB05 

2200 

LDFOO 

4777 

4777 

II-3 
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0157 

3445 

SBD45 

0160 

6653 

PJB53 

0161 

2044 

LDD44 

0162 

3445 

SBD45 

0163 

6003 

ZJF03 

0164 

6657 

PJB57 

0165 

0002 

ERR02 

0166 

7101 

JFI01 

0167 

0172 

WRITE 

0170 

1000 

0171 

7700 
ADDITIONAL  CELLS  USED 

HLTOO 

0044 

LIMIT 

(PUT  IN  LIMIT) 

0045 

ADDR 

7600 

0046 

CNTR  1 

0047 

CNTR  2 

0050 

BUFFER  STORAGE 

0051 

BUFFER  STORAGE 

0052 

BUFFER  STORAGE 

0053 

BUFFER  STORAGE 

0054 

BUFFER  STORAGE 

0055 

BUFFER  STORAGE 

0057 


NUMBER  OF  MEMORY  DUMPS 


II-4 
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MEMORY  DUMP  TO  M.T. 


H-5 


0172 

0173 
0174 

0175 
0176 

0177 
0200 
0201 
0202 

0203 
0204 

0205 
0206 

0207 
0210 
0211 
0212 
0213 
0214 
0215 
0216 
0217 
0220 
0221 
0222 
0223 
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5457 

A0D57 

2600 

LCF00 

7776 

4060 

STD60 

7500 

EXFOO 

1141 

REQ  STAT 

7600 

INAOO 

0202 

LPN02 

6004 

ZJF04 

5460 

A0D60 

6506 

NZB06 

0003 

EURO  3 

7500 

EXFOO 

2111 

2111 

7303 

0UT03 

7600 

TERM  OUT    (LI 

6102 

NZF02 

1000 

INIT  WRIT 

7500 

EXFOO 

1141 

REQ.  STAT 

7600 

INAOO 

0240 

LPN40 

6103 

NZF03 

7101 

JFI01 

0100 

00 

7700 

HLTOO 

DUMP  TO  PAPER  TAPE 


0163 

7500 

0164 

4104 

0165 

7303 

0166 

7600 

0167 

6102 

0170 

1000 

0171 

7700 

II-6 
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DAILY 

IDENTIFICATION 

USE  8  CELLS 

0030 
0031 

0070 
0064 

0032 
0033 
0034 

0032 

0030 

0060 

0035 
0036 

0072 
0062 

EXAMPLE 

1  day 

3 

J 

A  month 
N 


4  year 
0037  0045  CAR  RET 


11=7 
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DAILY  RECORD  TO  M.T. 


■SfpAKV 


; 

uPD 
31 

iot>n 

>/ 

LDF 

y 

^ 

Oloo 

, * 

Ifito 

; 

/ 

APZ>£ 

00 

J/ 

AD£>£ 

\/ 

LDD 

u 

S8D 
Aj>P£ 

n-8 
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DAILY  RECORD  TO  M„T 


0224 
0225 
0226 
0227 
0230 
0231 
0232 
0233 
0234 
0235 
0236 
0237 
0240 
0241 
0242 
0243 
0244 
0245 
0246 
0247 
0250 
0251 
0252 

0253 
0254 

0255 
0256 

0257 
0260 
0261 
0262 
0263 
0264 

0265 
0266 
0267 
0270 
0271 
0272 

0273 
0274 


2030 
4100 
1000 
2031 
4100 
1001 
2032 
4100 
1002 

2033 
4100 

1003 
2034 
4100 
1004 

2035 
4100 
1005 
2036 
4100 
1006 
2037 
4100 
1007 
2200 
7700 
4100 
0100 
2200 
1010 
4045 
0400 
4145 
2044 
3445 
6003 
5445 
6506 
7101 
0172 
7700 


LDD30 
STIOO 

1000 
LDD31 
STIOO 

1001 
LDD32 
STIOO 

1002 

LDD33 
STIOO 

1003 
LDD34 
STIOO 

1004 

LDD35 
STIOO 

1005 
LDD36 
STIOO 

1006 
LDD37 
STIOO 


>DDR) 


LDFOO 

7700 
STIOO 

0100 
LDFOO 

1010 
STD45 
LDNOO 
ST  1 45 
LDD44 

SBD45 

ZJF03 

A0D45  (ADDR) 

NZB06 

JFI01 

0172  (DUMP) 
NORMAL  STOP 


II-9 


REWIND,  READ  AND 


DAILY  RECORD  SEARCH 


Alma 


n.k  i.  InWlA      J.  J. 


REWIND,  READ  AND 


DAILY  RECORD  SEARCH 


» 


! 


• 


Y£"S 
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Start 

Rewind  0275 

0400 

LDN00 

0276 

4043 

STD43 

0277 

7500 

EXFOO 

0300 

1141 

REQ  STAT 

0301 

7600 

INAOO 

0302 

0202 

LPN02 

0303 

6504 

NZB04 

0304 

7500 

EXFOO 

0305 

1161 

REWIND 

0306 

7500 

EXFOO 

0307 

1141 

REQ  STAT 

0310 

7600 

INAOO 

0311 

0202 

LPN02 

0312 

6504 

N2B04 

Read   0313 

7700 

HLTOO 

Start  0314 

7500 

EXFOO 

0315 

1141 

REQ  STAT 

0316 

7600 

INAOC 

0317 

0202 

LPN02 

0320 

6504 

NZB04 

0321 

7500 

EXFOO 

0322 

2111 

EXF11  READ 

0323 

7203 

0UT03 

0324 

7600 

TERM  READ 

0325 

6102 

NZF02 

0326 

1000 

INIT  READ 

0327 

7500 

EXFOO 

0330 

1141 

REQ  STAT 

0331 

7600 

INAOO 

0332 

0240 

LPN40 

0333 

6002 

ZJF02 

0334 

0004 

ERR04  (END  OF 

0335 

2043 

LDD43   TAPE) 

0336 

0701 

SBN01 

0337 

6004 

ZJF04 

0340 

6525 

NZB25 

Search  0341 

0401 

LDN01 

Start  0342 

4043 

STD43  (SEARCH) 

0343 

2100 

LDIOO 

0344 

1000 

1000 

0345 

3430 

SBD30 

0346 

6532 

NZB32 

0347 

2100 

LDIOO 

0350 

1001 

1001 

0351 

3431 

SBD31 

0352 

6536 

NZB36 

0353 

2100 

LDIOO 

11-12 
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0013 
N2 
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0354 

1002 

1002 

0355 

3432 

SED32 

0356 

6542 

NZB42 

0357 

2100 

LDIOO 

0360 

1003 

1003 

0361 

3433 

SBD33 

0362 

6546 

NZB46 

0363 

2100 

LDIOO 

0364 

1004 

1004 

0365 

3434 

SBD34 

0366 

6552 

KZB52 

0367 

2100 

LDIOO 

0370 

1005 

1005 

0371 

3435 

SBD35 

0372 

6556 

NZB56 

0373 

2100 

LDIOO 

0374 

1006 

1006 

0375 

3436 

SBD36 

0376 

6114 

NZF14 

0377 

2100 

LDIOO 

0400 

1007 

1007 

0401 

3437 

SBD37 

0402 

6110 

N2F10 

0403 

2100 

LDIOO 

0404 

1010 

1010 

0405 

6105 

NZF05 

0406 

0400 

LDNOO 

0407 

4043 

STD43 

0410 

7101 

JFIOI 

0411 

0077 

0077 

0412 

7101 

JFIOI 

0413 

0314 

TEST  TAPE  22  BLOCKS 

DATA 

45io 

LIMIT 

CELL 

44  - 

7600 

J*\J 

ADDR 

CELL  45 

3  4051s 

CNTR  2 

CELL 

21 

Nl 

% 

CELL 

22 

«!■ 

N2 

CELL 

23 

CNTR  1 

CELL 

20 

0526 

7700 

HLTOO 

0527 

0400 

LDNOO 

0530 

4021 

STD21 

0531 
0532 

2200 
1000 

LDFOO 
1000 

11-13 
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0533 

4045 

STD45 

0534 

0413 

LDN13 

0535 

4023 

STD23 

0536 

0400 

LDN00 

0537 

4020 

STD20 

0540 

2200 

LDFOO 

0541 

4051 

4051 

0542 

4022 

STD22 

0543 

2021 

LDD21 

0544 

0726 

SFN26 

0545 

6417 

ZJBl? 

0546 

2023 

LDD23 

0547 

4145 

STI45 

0550 

5445 

A0D45 

0551 

2022 

LDD22 

0552 

4145 

ST  1 45 

0553 

5445 

A0D45 

0554 

5422 

A0D22 

0555 

5420 

A0D20 

0556 

2200 

LDFOO 

0557 

4777 

4777 

0560 

3445 

SBD45 

0561 

6207 

PJF07 

0562 

2044 

LDD44 

0563 

3445 

SBD45 

0564 

6104 

N2F04 

0565 

5421 

A0D21 

0566 

7101 

JFI01 

0567 

0172 

0172 

0570 

2020 

LDD20 

0571 

0777 

SBN77 

0572 

6434 

ZJB34 

0573 

6525 

NZB25 

11=14 
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WRITE  TEST  TAPE 


22  BLOCKS  DATA 


S'PAKT 


»o 


/ 


' 


XZL 


HA 


S8)J 

n\ 

■ 

iPD 

CW? 

i 

1 

6Bi> 
A  DP? 

dh 

4777 

/ 

\ 

AOJ> 

6A/A?^y 

L    i 

\ 

AOJ) 

4r- 

11-15 
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FORTRAN  PROGRAMS  PAGE 


ANALYSIS  PROGR'Lli 


SATELLITE  DISPLAY  PROGRAM 


III 


DIMENSION  NDATA (664 ), DATE (  O0),H(  17 28 ) , HONE ( 90C )  , ITITLE { 12)  , 
1HTWO(900),HTHRE(900) , T { SCO ) , B IN M  10 0  )  , E  I  N2 (  1 00 )  , B I N3 (  1 OC )  ,  X { 1 00  ) 
DO  5  J  =  1,14 
READ  5,DATE(J) 

5  F0RMAT(C16) 

CON (MASK  =  0C0CCC00  7777  777/B) 
ENA( NDATA) , INA( 1 ) . 
STA(  INI  T)  ,  INA(.C64)  . 
STA( I  TERM) ,LNI(OJ. 

6  CALL  READ(3, INITtlTERM, 1) 
DO  70  I=3,.; 

SLS1  (7) 

7  ADATE  =  DATE( I ) 
BDATE  =  DATEU+1  ) 

LD A (ADATE) ,SUR(NDATA  +1),AJP1(6), 

LD A (BDATE) , SUB (NDATA  +2  ) , AJP1 ( 6 ) . 

NBLGC  =  0 
30  CALL  READto, INIT, ITERM,  1) 

N2RC  =  0 

DO  81  J=l ,864 

IF(NDATAU))  82,82,81 
82  NZRO  =  NZRG+1 
81  CONTINUE 

ND  =  NDATA(3) 

LDA(ND) ,AJP0( 10 ). 

K  =  C 

DO  iJ  0  J  =  l  ,864 

K  =  K+1 

INDATA  =  NDATA(J) 

LDA( INDATA) , LOG (MASK) ,ARS(2U) ,STL(INLIST). 

H(K)=INLIST 

K  =  K  +  1 

LDA( INDATA) , LOG ( MASK ), STL ( 1NLIST) . 

SLS1 (40). 
40  H(K)=INLIST 

AVt=G. 

DO  50  J  =1, 1728 

H(J)  =  H(J)/3. 

H( J)=1024./H( J) *2348400.0 

50  AVE  =  H(J)  +  AVE 

PRINT  26,  AVE 
C  ZAPP  FILTER 

NERR  =   0 

AVE  =  0. 

HBAR  =  50950.0 

DO  51  J=l,1728 

DEV  =  ABSF(H( JJ-HBAR) 

IFIDEV-200.)  51,52,52 

JCGUNT  =  J 

52  NERR  =  NERR  +  1 
IFU-3)  53,53,54 

53  H(J)  =  H(J+5) 
GO  TG  51 

54  HtJ)  =  H(J-3) 

51  AVE  =  AVE  +  H(J) 
AVE  =  AVE/1728. 
PRINT  26,  AVE 

PRINT  97, NERR, NZRO, JCOUNT 

NBLOC  =  NBLOC  +  1 

WRITE  TAPE  2,  H 

GO  TO  50 
10  REWIND  2 
C  THIS  SECTION  COMPUTES  H  VS  TIME 
C  18  PT  AVE  N  IS  TOTAL  NUMBER  GF  POINTS  PER  DAY 

TIME  =  -  6.5/704. 

DO  60  L=l, NBLOC 

READ  TAPE  2.H 

DO  62  J=l,32 

Hi  BAR  =  0. 

H2BAR  =  0. 

H3BAR  =  0. 

DO  61  K=l  ,18 

M  =  5*K  -  2  +  (J  -  l)«5l4         III-l 

H1BAR  =  H(M)  +  HlBAR 

H2BAR  =  H(M+1 )  +  H2BAR 

61  H3BAR  =  H(M+2)  +  H3BAR 
N  =  J  +  (L-l )*3  2 
HONE(N)  =  Hi  BAR/18. 
HTWC(N)  =  H2BAR/1-8. 
HTHRE(N)  =  H3BAR/18. 
TIME  =  TIME  +  6.5/704. 

62  T(N)  =  TIME 

.   .    L  r\       C  f>  "k.'-T  T  Mi  t  C  - . „-_ 


REMIND    2 
PR  1  N'T    2  7,  N 
PRINT    27, M 
NPTS=N 
MEAN    VALUE    OF    DIFFERENCES 
JERR    =    0 
I  ERR    =    0  APPENDIX  in 

NERR  =  0  

DELI  =  0. 

DEL2  =  0. 

DEL3  =  0. 

DO  61;  J  =  1  ,NBLCC 

READ  TAPE  2,H 

DO  64  L  =  1,576 

M  =  3*L-2 

DEL  =  H(M  +  1 )  -  H(M) 

IFtARSF(DEL)-  UO. )  65,66,66 

65  DELI  =  DELI  +  DEL 
GO  TO  67 

66  NERR  =  NERR  +  1 

67  DEL  =  H(M  +  1  )  -H(M  +  2  ) 
IFtABSF(OEL)  -  40.)  68,69,69 

68  DEL2  =  DEL2  +  DEL 
GO  TO  71 

69  I  ERR  =  IERR  +  1 

71  DEL  =  H(M)  -  H(M  +  2) 
IF(ABSF(DEL)  -  40.)  72,73,75 

72  DEL3  =  DELi  +  DEL 
GO  TO  64 

73  JERR  =  JERR  +  1 
64  CONTINUE 

REWIND  2 

PRINT  97, NLRR, IERR, JERR 
97  F0RMAK3I20) 

ABLCC  =  NBLOC 

DELI  =  DELI/ (576. *AELOC) 

DEL2  =  DEL2/ (576.*A8LCC  ) 

DEL3  =  DEL5/ (576.»ABL0C) 

PRINT  49,DEL1,DEL2,DEL3 
DENSITY  FUNCTION 

STUEV1=0. 

ST0EV2=0. 

STDEV3=0. 

NERR  =  0 

IERR  =  0 

JERR  =  0 

XL=-25. 

DO  4  5  K=l,100 

BIN1 (K)=C. 

BIN2(K)=0. 

EIN3(K)=0. 

X(K)=XL 
45  XL  =  XL+0.5 

DO  74  J  =  l , NBLOC 

READ  TAPE  2»H 

DO  74  L  =  1,576 

M  =  3*L-2 

STOEVl=  (H(M+1 )-H(M)-DELl)*«2  +STDEV1 

STDEV2=  (H(M+1 )-H(M+2)-CEL  )**2  +  ST0EV2 

STDEV3=  (H(M)-H(M+2)-0fcL3)«*2  +STDEV3 

DEL  =  H(M+1 )  -  H(M) 

IF (ABSF (DEL-DEL  I )-25. )75,7i,76 

76  NERR  =  NERR  +  1 

75  DEL  =  H(N+1 )  -H(M+2) 

IF( ABSF(DEL-CEL2)-25. ) 7 7, 77, 78 
78  IFRR  =  IERR  +  1 

77  DEL  =  H(M)  -  HlK+2) 

IF( ABSF( DEL-DEL  3) -25. ) 2 9, 29,80 
80  JERR  =  JERR  +  1 
29  DO  31  K=l ,3 

IF(K-2)32,33,34 

32  D£L  =  H(.Y+1  )-H(M)-D£Ll 
GO  TO  55 

33  DEl=H(M+7)-H(M*2)-DEL2 
GO  TO  55 

34  DEL  =  H(IM)-H(M  +  2)-DEL3 

35  XL=  -25. 

DO  36  N=l ,100 

IF(DEL-XL)  37,37,38  III-2 

38  IF(DEL-(XL+0.5)  )39,39, 57 

39  IF(K-2)  42,43,44 

42  PIN! (N)=B INI (N)+l . 
GO  TO  31 

43  BIN2(N)  =  BIN2(N1+1 . 
GO  TO  31 


c 


NORMALIZED       DENSITY 
AMAX=BIN1 ( 1 ) 


APPENDIX  III 


206 


*4=*  +  //**********S  // 

CONTINUN 

PRINT    49, AMAX, BMAX.CMAX 

DO  206  J=l, TOO 

TINT ( J)=BIN1 ( J)/AMAX 

>IN2( J)=BIN2( J)/PMAX 

UN3(J)*BIN3{J)/CMAX 

)0  46  K  =  l,  12 


>3 

M 1 00 , X  ,  B I N 1 , 0 , C , L  A  B A , I T I TL  E , 0 , 0 , 0 , 4 , 0 , 2 , 8 , 8 , 0 , LA  ST ) 
..M100,X,BIN2,0,C,LABB,  ITITLE,0,0,0,4,0,2,8,8,C,LAST) 
CALL  DRAW { 1 00, X, 8 1N3, 0, C, LA 6C, I  TITLE, 0,0, 0,4, 0,2, 8, 8,0, LAST) 
506  N=NPTS 

HONEAV=0. 
HTWCAV=0. 
HTHREAV=C. 
DO  48  L=1,N 
HONEAV=HONEAV+HCNE( L) 
HTWCAV=HTWi,AV+HTWC(  L  ) 
HTHREAV=HTHREAV+HTHRE(L) 
SV04  =  N 

HONEAV  =  H0NEAV/SV04 
HTWCAV  =  HTWCAV/SV04 
HTHREAV  =  HTHREAV/SV04 
PRINT  49, HONEAV, HTWOAV, HTHREAV 
FORKAK 3E20.7) 
PRINT  27,NBLGC 
FORMATt 17) 

HONEAV*  (HONEAV  +HTWOAV  +HTHREAVJ/3. 
PRINT  26, HONEAV 
26  FOR?'AT(  E20.7) 

WRITE  TAPE  2,  HONE 

WRITE  TAPE  2,  HTWO 

WRITE  TAPE  2,  HTHRE 

REWIND  2 

DO  90  L  =  1  ,  N 

HONE(L)  =  HONE(L)  -  HONEAV 

uti.jpi  —    uTuni  i  i     _    HONEAV 

-  HCNEAV 


48 


49 
27 


HONE(L)  =  HONE(L)  -  HONEAV 
HTWC(L)  =  HTWO(L)  -  HCNEAV 
90  HTHRE(L)  =  HTHRE(L)  -  HCNE 
DO  1  L=l,12 
1  ITITLE(L)  =  8H 


III-3 


LABC=4H  H3 

CALL  DR AW (N, T, HCNE, 1,0, LAB A, ITITLE,  1.0, 10. 0,5,  I, 2, 2, 8, 10,0, LA  ST) 
CALL  DRAW (N, T, HTWO, 2^0, LABC,  I  TITLE,  1  . 0 , 1 0  .  C  ,  5 , 1  2  2  8 , 1 0, 0 , LA S T 
CALL  DRAW(N,T,HTHRE,3,0,LABC,ITITLE, 1 . C , 1 0. ,5, 1  J  2 ,2 , 3, 1 c! 0 , LAST 
ANSFORM  LOW  PASS  FILTER  N  ECUAL  20  TITUS                         ' 
at-t  i,  «.?  ),  111  u/inr, 


C    L  I  Krtli  Jl 

at-  i,  »g  i.ii,  u/mn 


IFU-2C)    302,303,303 
502    M=J 

GO    TO   30U 
3  03    M=20 
304    DO    301     K=1,M  APPENDIX  III 

L=    K-l  

BL=L 

ASUM 

BSUM 

CSUK 

OS  UN 


NLAGS=75 

MPTS  =    N-NLAGS 

PTS  =    MPTS 

XL  =  0. 

DO    9    J=1,100 

BINl ( J)=0. 

BIN2<J)«0. 

BIN3( J)=0. 

X( J)=XL 

XL    =    XL+1 


READ    TAPE    2, HONE 

READ    TAPE    2,HTW0  TTT-i, 

-,  HTHRE  J"U"   4 


READ    TAPE 
REWIND    2 

do  n  j=i 

H 

H 

H 

DO 

DO 

L=K-1 


i  i<->m  r-  /  \/   \     -  i  i  f*>  I  • 


EliMJlJ)=  blN^(J)  +  HTWO(K)*HTWO(L+J ) 
14  BIN.-5C  J)=BIN3  {  J)  +  HTHRt(K)*HTHRL(L+J  ) 

BINl ( J)=BIN1 ( J)  /PTS 

BIN2(J)=  8IN2(J)/PTS 
12  R I r i 3 ( J ) =  BIN3U)/PTS 
SAVE  AUTOCORRELATION   FUNCTION 

00    105    J=l,100 

H(J)=C1N1U)  APPENDIX  III 

H( J+100)=BIN2( J)  

105  H{ J+200)=BIN3( J) 
00    100    J=  1 , 1 2 

100  ITITL£(J)=-H 
ITITLE(  1  )=8hANt)ERSCN 
ITITLE(2)=oH  BOX  2^3 
ITITLE(7)=8HAUTCC0RR 
ITITLE(8)=FHELAT  ION 
ITITLE(9)=3HFUNCTI0N 
ITITLE( 10)=0H  Y  IN  PR 
ITITLE( 11 )=8H00UCTS  X 
ITITLE( 12)=8H  IN  LAGS 
LABA=4H  HI 

LABB=4H  H2 
LABC=4H  H3 

CALL  DRAW(NLAGS,X,BIN1,  1,0,LABA, I  TITLE, 0,0,0,0,0, 0,8,8,0, LAST) 
CALL  DR AW( NL AG S,X, 81 N2, 2, 0,L ABB, I  TITLE, 0,0,0,0,0, 0,8, 8,0, LAST) 
CALL  OR  AW ( NL AG S, X, BI N3, 3,0, LA BC, I  TITLE, 0,0, 0,0, 0,0, 8, 8,0, LAST) 
CROSS  CORRELATION  FUNCTION  75  LAGS 
DO  101  J= 1,100 
BIN1(J)=0. 
BIN2( J)=0. 

101  PIN3(J)=0. 

DO  102  J=1,NLAGS 

DO  103  K=1,MPTS 

L=  K-l 

BIMl ( J)=PIN1 ( J)+  HTWO(K)«HnNE(L+J ) 

B  I  s\  2  {  J)=8IN2  (  J)  +HTWO(K)*HTHRE{L+J) 

103  BIN3( J)=BIN3( J)+  HONE ( K )*HTHRE ( L+J) 
BIN1 ( J)=BIN1 { J)/PTS 
8IN2(J)=  BIN2(J)/PTS 

102  BIN5(J)=  BIN3(J)/PTS 
DO  10M  J=l,12 

104  ITITLE(J)=8H 
ITITLE( 1 )=8HANDERS0N 
ITITLE(2)  =8H  BOX  263 
ITITLE(7)  =6HCR0SC0RR 
ITITLEC8)  =  8HEL/UI0N 
ITITLE(9)  =0HFUNCTION 
ITITLE(  10)=6H  Y  IN  PR 
ITITLEt 11 )=8H0CUCTS  X 
ITITLE( 12)=8F  IN  LAGS 
LA6A  =  UH  H21 
LABB  =  UH  H23 
LABC  =  UH  H13 

CALL  DRAW(NLAGS,X,BIN1,  1, 0, LA BA, I  TITLE, 0,0, 0,0, 0,0, 8, 8,0, LAST) 
CALL  DR AW ( NL AGS, X, B I N2, 2,0, LABB, I  TITLE, 0,0, 0,0, 0,0, 8, 8, 0,LAST) 
CALL  DRAW(NLAGS,X,BIN3, 3,0, LABC, I TITLE, 0,0, 0,0, 0,0, 8, 8,0, LAST) 

COHERENCY   COMPUTATION  75  LAGS 
DO  106  J  -    1 ,NLAGS 

BIN1 (J}s(BINl( J)»*2)/(H(J+100)«H(J) ) 
IF(BINKJ))  107,108,108 

107  BIN! { J)=-BIN1 (J) 

108  BIN! ( J)=  SCRTFI BIN!  ( J)  ) 
BlN2(J)=(BIN2(J)*«2)/(H(J+i00)»H(J+200)) 
IF(BIN2(J))  109,110,110 

109  BIN2( J)=-8IN2( J) 

110  BIN2(J)=  SORTF( EIN2( J)  ) 
BIN3(J)=(BIN3{J)**2)/(H(J}*H{J+200> ) 
IF(BIN3(J))  111,106,106 

111  BIN3( J)=-BIN3{ J) 

106  BIN3(J)=  SLRTFt BIN3( J) ) 
ITITLE(7)=8HC0HERENC 
ITITLE(8)=8HE 
ITITLE( 10)=8H  Y  IN  CO 
ITITLEt 1 1 )=RHH  X  IN  L 
ITITLEt 12)=8HAGS 

CALL  CRAW(NLAGS,X,BIN1,1,0,LABA,ITITLE,0,.2,0,0,0,0,8,8,0,LAST) 
CALL  DRAW(NLAGS,X,BIN2,2,0,LABB, I TI TLE , 0, . 2 , 0, C , 0 , 0 , 8 , 8 , 0 , LAS T ) 
CALL  DRAW(NLAGS,X,BIN3, 3,0,LABC,ITITLE,0,.2,0,0,0,0,8,8,C,LAST) 

70  CONTINUE 

REMIND  3  III-5 

STOP 

END 

MACHINE  READ(1ARG,2ARG, 3ARG,4ARG) 


*  TYPEWRITER  MESSAGE  CCOE  TABLE 

CON(T6  =U504CU122050220UP, 

*  TB  =121203]2COOOOCOCP, 

*  no  =  1530121U01250U2CE, 

CON  (  Tl  2  =  4 50U 04  1220  ^0220*4  0, 

*  Tiu=iui 1 20OOOO0OO0OCP, 

*  TU0=450404013015200M[', 

*  T42=C41406221416300 IE, 


ArxnnvjJxA    XXX 


,  121 401 20 B, 
T4  ;  =340  6  140 1040 6C3C 18, 
T43=202  2  420000CCC0C0B) 


TAPE  UNIT  ASSIGNMENT  TABLE 

ASSIGNMENTS  FOR  CHANNEL  3/4  LIBRARY  TAPE 


CON(KO  =0  , 
«  K3  =3320306, 
*     K6  =5520206, 


Kl  =3220  OB, 
K4  =  3220*108, 
K7  =5520308, 


K2  =3320206, 
K5  =552010  6, 
K8  =55^GU08) 


*  ASSORTED  CGNSTANTS 

CONCRU  =73737373737373738, 


Rl  6  =  007 OCCCCCOOOOOOOB) 


»  TAPE  RtADING  ROUTINE 


IX 

1ARG 
2ARG 
3ARG 
UARG 


2X 

3X 

4X 
5X 
6X 


7X 
BX 
9X 


10X 

1  IX 
12X 

13X 
14X 
15X 
17X 


J  oX 


19X 


SLJ( 
ZRi.'{ 
ZRO( 
ZRiM 
ZRO( 
L0A7 
SIU1 
ENA( 
SLJ( 
INA{ 
LIL1 
ADD7 
SCL( 
SAL( 
INA( 
INA  ( 
INA{ 
INA( 
LRSt 
LLS( 
SAL( 
LLS( 
STA{ 
SSU( 
LDA7 
LDA( 
STA{ 
LDA7 
SAU( 
ENI2 
EXF( 
EXF( 
EXF( 
SU«{ 
IN  I  1 
INI  1 
EXF7 
EXF7 
EXF7 
QJPl 
SLJ( 
I  JP2 
IJP2 
AJP{ 
EXF  { 
LILT 
ENA( 
ENA( 
ENA( 
SLJ4 
AJP  ( 
SLJ  i 
SLJU 
AJP( 
SLJ( 
ENI  1 
SLJ( 


N) 

0) 

0) 

0) 

C) 

(1ARG) 

(19X) 

TuO) 

P) 
116) 
(KO) 
(4ARG) 
7776) 
5X) 
3) 
2) 
1) 
1  ) 
18) 
3) 
3X) 
39) 
KO) 
KO) 

(2ARG) 
6X) 
6X) 

(3ARG) 
6X) 
(4) 
N) 
N) 
N) 
R4) 
(-1  ) 
(1  ) 
(N) 
(N) 
(N) 
(  19X) 
19X) 
(L+2) 
(L+  1  ) 
3X) 
N) 
(N) 
T6) 
T9) 
T12) 
(  G) 
19X) 

nx) 

(  G) 

3X) 
11X) 
(N) 
IX) 


SLJ( 
ZRC( 
ZRO( 
ZRO( 
ZRO( 
INA( 
AJP3 
SLJU 
ZRO( 
STA( 
LDAl 
SAU( 
SAL( 
SAL( 
SAU( 
SAU( 
SAU( 
SAU( 
ENA( 
SAL( 
SAU( 
LDQ( 
LGA{ 
STA( 
SAU( 
SSU( 
ENH 
SAL( 
SIL2 
LIL1 
EXF7 
EXF7 
LDA7 
LLS( 
L0A1 
SUB( 
SLJ( 
SLJ( 
SLJ( 
AJP{ 
ZRC( 
SLJ( 
SLJ{ 
QJP( 
EXF7 
SLJ( 
SLJ{ 
SLJ( 
SLJ( 
ZRO( 
AJP3 
ZRO( 
ZRfH 
AJP3 
ZRO( 
EN  I  2 
ZRO( 


) 


L+5) 

0) 

0) 

C) 

0) 

-1  IB) 

(L+3) 

(  G) 

0) 

KO  ) 

(KC) 

4X) 

4X  )  . 

11X) 

5X) 

9X) 

nx 

7X) 

0) 

2X  ) 

12X) 

RU) 

5X) 

5X) 

5X  ) 

KO  ) 

0) 

N) 

(  19X) 

(N) 

(N  ) 

<N) 

(2ARG) 

US) 

(I) 

RU  J 

13X) 

L  +  M 

10X) 

3X) 

0) 

1UX) 

13x) 

3X  ) 

<N) 

4X) 

17X) 

17X) 

18  X) 

0) 

(3X) 

0) 

0) 

( 3  X  J 

0) 

(N) 

0) 


III-6 


EXIT/ENTRY 

TAPE  UNIT  NUMBER  ARG. 

INITIAL  ADDRESS 

TERMINAL  ADDRESS 

MODE,  1  =  BINARY,  2  =  BCO 

CHECK  FOR  TU  ASSIGNMENT 
GO  PRINT  ERROR 


SET  TU  CONTROL  CODES 
INITIALIZE 


SELECT  UNIT 
ACTIVATE  CHANNEL 


SENS 
SENS 
SENS 
SENS 
EXIT 
SENS 
SENS 
SENS 
BACK 
RETU 
SET 


E  ENDFILE 
E  PARITY  ERROR 
E  LENGTH  ERROR 
E  FOR  BAD  RECORD 

E  FIVE  PARITY  ERRORS 

E  FIVE  LENGTH  ERRORS 

E  FOR  BAD  RECORD 

SPACE 

RN  TO  REREAD 

TO  INDICATE  TROUBLE 


TROUBLE  EXIT 

SENSE  ACTION  TO  TAKE 

TROUBLE  EXIT 

SENSE  ACTION  TO  TAKE 

RE-STORE  INDEXES 


■  ■  —--■.-..,  ■  —  ■>, 


PROGRAM    PROMAG  APPENDx^  III 

DIMENSION    NDATA ( 864 ), DATE ( 100),H( 1728 ), HONE ( 900 ) , ITITLE( 12) , 
1HTW0(900) ,HTHRE(900) ,T(900) 
DO    5    J    =     1  ,  14 
READ    5,DATE(J) 

5  F0RMAK016) 

CONtMASK  =  0000000077777777B) 
ENA(NDATA) , INA( 1 ). 
STA( INIT)  , INA(864). 
STA( ITERM),ENI (0) . 

6  CALL  READ(3, INIT, ITERM,  1) 
DO  70  I  =  1,11,2 

ADATE  =  DATEU  ) 

BDATE  =  DATE( 1  +  1  ) 

LDA( ADATE) ,SUB(NDATA  +1),AJP1(6), 

LDA( BDATE) ,SUB(NDATA  +2),AJP1 {6). 

NBLOC  =  0 
30  CALL  READ(3, INIT, ITERM,  1) 

ND  =  NDATA(3) 

LDA(ND) ,AJP0(10). 

K  =  0 

DO  40  J=l,864 

K=K+1 

INDATA  =  NDATA(J) 

LDA( INDATA) , LDQ ( MASK ) , ARS { 24 ) , STL { I NLI ST ) . 

H(K)=INLIST 

K=K+1 

LDA( INDATA) , LDQ (MASK) , STL { INL 1ST ) . 

SLS1 (40). 
40  H(K)=INLIST 

AVE=0. 

DO  50  J  =1,1728 

H( J)=( 1024./H( J) )*23484C0.0 

50  AVE=H(J)+AVE 
AVE=AVE/1728. 
PRINT  26,  AVE 

ZAPP  FILTER 

NERR  =   0 

AVE  =  0. 

HBAR  =  50950.0 

DO  51  J=l,1728 

DEV  =  ABSF(H{ JJ-HBAR) 

IF(DEV-200.)     51,52,52 

52  NERR  =  NERR  +  1 
IFU-3)    53,53,54 

53  H( J)  =  H( J+3) 
GO  TO  51 

54  H(J)  =  H(J-3) 

51  AVE  =  AVE  +  H(J) 
AVE  =  AVE/1728. 
PRINT  26,  AVE 
PRINT  27,  NERR 
NBLOC  =  NBLOC  +  1 
WRITE  TAPE  2,  H 
GO  TO  30 

10  REWIND  2 

NPTS    =    (  1728«NBLOC)/2700 

PRINT  27, NPTS 

READ  TAPE  2,  H 

K  =  1 

TIME  =  -6.5/900. 

M=l 
THIS  SECTION  COMPUTES  H  VS  TIME 

DO  60  L=l,900 

H0NE(L)=(H(K)+H(K+3)+H( K+6 > )/3. 

HTWC(L)  =  (H(K  +  1  )+H(K  +  4)  +  }-:(K  +  7)  )/3. 

LCOUNT=L 

HTHRE(L)=(H(K+2)+H(K+5)+H(K+8) ) /3. 

TIME  =  TIME  +  6.5/900. 

T(L)  =  TIME 

K  =  K  +  3*NPTS 

MAX  =  K+8 

IF(MAX  -  1728)  60,60,80 
80  K=l 

IF(M-NBLOC)  85,95,95  TTT  7 

85  READ  TAPE  2,H  J"LX  ' 

M=M+1 
60  CONTINUE 
95  ?RIKT    97-,  LCOUNT 

PF INT  97  M 

AT    ^    4,*3.1416/100. 

REWIND    2 
.  WRITE  TAPE  2,  HONE 

WRITE  TAPE  2,  HTWO 

i.iottc     t  a  r>  r     ^  uTimc , , __, _, 


riuiN  l_  m  v  —  U  • 

HTWCAV=0. 

HTHREAV=0. 

DO    U8    L=l,900 

HONEAV=HONEAV+HONE(L)  APPENDIX  III 

HTWCAV=HTWOAV+HTWO(L)  

U8  HTHREAV  =  HTHREAV  +  HTHRE (L  ) 

H0NEAV=H0NEAV/900. 

HTW0AV=HTWOAV/900. 

HTHREAV=HTHREAV/900. 

PRINT    49, HONEAV, HTWOAV.HTHREAV 
U9    F0RMAK3E20.7) 

PRINT  27,NBL0C 
27  FORMATt 17) 

HONEAV=  (HONEAV  +HTWOAV  +HTHREAV)/3. 

PRINT  26, HONEAV 
26  FORMAT(E20.7) 

DO  90  L=l ,900 

HONE(L)  =  HONE(L)  -  HONEAV 

HTWC(L)  =  HTWO(L)  -  HONEAV 
90  HTHRE(L)  =  HTHRE(L)  -  HCNEAV 

DO  1  L=l,12 
1  ITITLE(L)  =  8H 

ITITLE(  1  )=8HANDERS0N 

ITITLE(2)=8H  BOX  263 

ITITLE(7)=8H  EARTHS 

ITITLE(8)=8HMAGNCTIC 

ITITLE(9)=8H  FIELD  V 

ITITLEt 10)=8HS  TIME 

ITITLE( 1 1 )=8HT  IN  HRS 

ITITLEt 12)=8H  H  GAMMA 

LABA  =  liH  HI 

LABB=4H  H2 

LABC=UH  H3 

CALL  DRAW(900,T,H0NE, 1,0,LABA,ITITLE,1.0,10.0,5,1,2,2,8,10,0,LAST) 

CALL  DRAW(900,T,HTWO,2,0,LABB,ITITLE,1.0,10.0,5,1,2,2,8,10,0,LAST) 

CALL  DRAW(900,T,HTHRE,3,0,LABC,  ITIFLE, 1.0, 10. ,5, 1,2, 2, 8, 10,0, LA  ST) 

CALL  S ATGR AF ( 0, 900, T, HONE, 0,0, 0,0) 

CALL  S ATGR AF { 0, 900, T, HTWO, 0,0, 0,0) 

CALL  SATGRAF(0,900,T,HTHRE,0,0,0,0) 
C  THIS  SECTION  COMPUTES  SIMPLE  DIFF 

READ  TAPE  2,  HONE 

READ  TAPE  2,  HTWO 

READ  TAPE  2,  HTHRE 

REWIND  2 

DO  15  L=l ,900 

H(L)=HONE(L) 

HONE(L)=H(L)-HTWO(L) 

STCRE=HTWO(L) 

HTWO(L)=H(L)-HTHRE(L) 

15  HTHRE(L)=STORE-HTHRE(L) 
DO  16  L=1,12 

16  IT1TLE(L)=8H 
ITITLEt 1 )=8HANDERS0N 
ITITLE(2)=8H  BOX  263 
ITITLE(7)=3HFIRST  DI 
ITITLE(8)=8HFFERENCE 
ITITLE(9)=8H  IN  MAG 
ITITLE( 10)=8HFIELD 
ITITLE( 1 1 )=8HT  IN  HRS 
ITITLE( 12)=8H  H  GAMMA 
LABA  =  l+HDELl 
LABB=UHDEL2 
LABC=UHDEL3 

CALL  DRAW(90  0,T,HONE, 1,0,LABA,ITITLE,1.0,2.00,5,1,2,2,8,10,0,LAST) 
CALL  DRAW(900fTfHTW0t  2, CfL ABB, ITITLEt 1.0.2.00,5, 1,2, 2f 8 ,10,0 f LAST) 
CALL  DRAW(900,T,HTHRE,3,0,LABC, ITITLE, 1.0,2.0,5, 1,2, 2, 8, 10,0, LAST) 
CALL  SATGRAF(0,900,T,H0NE,0,0,0,0) 
CALL  SATGRAF { 0, 900, T, HTWO, 0,0, 0,0) 
CALL  S ATGR AF ( 0, 900, T, HTHRE, 0,0, 0,0) 
70  CONTINUE 
REWIND  3 
STOP 
END 
MACHINE  READM  ARG,  2ARG,  3ARG  ,  UARG  ) 

*  NOTE   G  =  ADDRESS  OF  ENTRY  TO  TYPEWRITER  SENSE  ROUTINE 

*  P  =  ADDRESS  OF  PROGRAM  CONTROL  ROUTINE  JUMP   ENTRY 

*  THIS  ROUTINE  IS  FOR  PROGRAM  CALL 

L0C(G  =  31,   P=10)  II]C'3 

* 

»  TYPEWRITER  MESSAGE  CODE  TABLE 


*  TAPE  UNIT  ASSIGNMENT  TABLE 

*  ASSIGNMENTS  FOR  CHANNEL 


T1 3=200! 

T  1  5  =  4  50  4 u 4.1  1  \d.  itui^us, 
T41 =  340 6 140  104 060301 B. 
T4 3=202  24  200000000006) 


3/4  LIBRARY  TAPE 


CON(KO  =0  , 
«  K3  =3320300, 
•     K6  =5520206, 


K1  =3320106, 
K4  =3320406, 
K7  *5520iOB, 


K2  =3320200, 
K5  =5520106, 
K8  =5520400) 


*  ASSORTEO  CONSTANTS 

C0N{R4  =73737373737373738, 


R 16=007  00000000000006) 


«  TAPE  READING  ROUTINE 


IX 

1ARG 

2ARG 

3ARG 

4ARG 


2X 

3X 
4X 
5X 
6X 


7X 

8X 
9X 


10X 

11X 
12X 
13X 
14X 
15X 
17X 


18X 
19X 


29 


SLJ(N) 
ZRO(O) 
ZRO(O) 
ZRO(O) 
ZRO(O) 
LDA7(1 A 
SI Ul (19 
ENA(T40 
SLJ(  P) 
INA(  118 
LIL1 (KO 
ADD7(4A 
SCL(777 
SAL(5X) 
INA(3) 
INA(2) 
INA(1 ) 
INA( 1 ) 
LRSt 18) 
LLS(3) 
SAL(3X) 
LLS(39) 
STA(KO) 
SSU(KO) 
LDA7J2A 
LDA(6X) 
STA(6X) 
LDA7(3A 
SAU(6X) 
ENI2(4) 
EXF(N) 
EXF(N) 
EXF(N) 
SUB(R4) 
INIK-1 
INI1 (1 ) 
EXF7(N) 
EXF7CN) 
EXF7(N) 
QJP1 (19 
SLJC19X 
IJP2(L+ 
IJP2(L+ 
AJP(3X) 
EXF(N) 
LIL1 (N) 
ENA(T6) 
ENA(T9) 
ENA(T12 
SLJ4(  G 
AJP( 19X 
SLJ(11X 
SLJ4(  G 
AJP(3X) 
SLJ(  11X 
ENI1 (N) 
SLJ(IX) 
END 
END 
.51430 
50930 


RG) 

X) 

) 

) 
) 

RG) 
6) 


RG) 
RG) 


X) 
) 

2) 
1) 


83E+05 
95E+05 


SLJCL+ 

ZRO(O) 

ZRO(O) 

ZRO(O) 

ZRO(O) 

INA(-1 

AJP3(L 

SLJ4( 

ZRO(O) 

STA(KO 

LDA1 (K 

SAU(4X 

SAL( 4X 

SAL( 1 1 

SAU(8X 

SAU(9X 

SAU( 11 

SAU(7X 

ENA{0) 

SAL(2X 

SAU( 12 

LDCMRl 

LDA(5X 

STA(5X 

SAU(5X 

SSU(KO 

ENKO) 

SAL(N) 

SIL2( 1 

LILKN 

EXF7(N 

EXF7(N 

LDA7(2 

LLS(48 

LDA1 (Z 

SU8(R4 

SLJ( 15 

SLJ(L+ 

SLJ(  10 

AJP(3X 

ZRO(O) 

SLJ( 14 

SLJ(  13 

QJP(3X 

EXF71N 

SLJ(4X 

SLJ(  17 

SLJ(  17 

SLJ(  18 

ZRO(O) 

AJP3(i 

ZRO(O) 

ZRO(O) 

AJP3(3 

ZRO(O) 

ENI2(N 

ZRO(O) 


5) 


18) 
+  3) 
G) 

) 

C) 

) 

) 

X) 

) 
) 

X) 

) 
) 

X) 
6) 

) 
) 
) 
) 


9X) 

) 

) 

) 

ARG) 

) 

) 

) 

X) 

U) 

X) 

) 

X) 
X) 

) 
) 

) 

X) 
X) 
X) 

X) 


X) 

) 


EXIT/ENTRY 
TAPE  UNIT  NUM8ER 
INITIAL  ADDRESS 
TERMINAL  ADDRESS 
MODE,  1  =  BINARY 


ARG. 


2  =  DCD 


CHECK  FOR  TU  ASSIGNMENT 
GO  PRINT  ERROR 


SET  TU  CONTROL  CODES 
INITIALIZE 


SELECT  UNIT 
ACTIVATE  CHANNEL 


SENSE 

SENSE 

SENSE 

SENSE 

EXIT 

SENSE 

SENSE 

SENSE 

8ACKSPACE 

RETURN  TO  REREAD 

SET  TO  INDICATE  TROUBLE 


ENDFILE 
PARITY  ERROR 
LENGTH  ERROR 
FOR  BAD  RECORD 

FIVE  PARITY  ERRORS 
FIVE  LENGTH  ERRORS 
FOR  BAD  RECORD 


TROUBLE  EXIT 

SENSE  ACTION  TO  TAKE 

TR0U6LE  EXIT 

SENSE  ACTION  TO  TAKE 

RE-STORE  INDEXES 

III-9 


.5097462E+05 

.cnoOA  ~3  c  r-  .  r\c 


APPENDIX  IV 


CIRCUIT  DIAGRAMS 
POLARIZING  AND  DAMPING  NETWORK 
MAGNETOMETER  PREAMP 
MAGNETOMETER  AMPLIFIER 
RADAR  GATE 

CLOCK  GATE  AND  GATE  CONTROL 
STORAGE  RESET  AND  GATE  DRIVER 
TIMER 

COUNTER  START  AND  RESET 
SIGNAL  GATE  AND  GATE  CONTROL 
CLOCK  FREQUENCY  MULTIPLIER 
COUNTER  FLIP-FLOP  AND  GATE 
MULTIPLEX  PULSE  GENERATOR 
MULTIPLEX  BLOCK  DIAGRAM 
ANALOG  LADDER 
OUTPUT  AMPLIFIERS 

PLUS  AND  MINUS  12  VOLT  POWER  SUPPLY 
PLUS  AND  MINUS     h  VOLT  POWER  SUPPLY 
PLUS  28  VOLT  OSCILLATOR  AND  TIMER  POWER  SUPPLY 
INTERFACE  DRIVER  AND  AMPLIFIER 
INTERFACE  AND/OR  GATES 
INTERFACE  FLIP-FLOP  AND  INVERTER 
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